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Abstract. Based on the free precession model of the isolated pulsar PSR B1828-11, 


Link & Epstein ( 2001 ) showed that the observed pulse durations require the radio beam 
to have a non-standard shape: the beam duration is larger for beam sweeps farthest 
from the dipole axis. In their analysis they assumed that the actual precession period is 
~ 500 d. Recent theoretical studies suggested that the actual precession period might be 
~ 1000 d as seen in observations (Rezania 2002, Wasserman 2002|). In this paper, in a 


good agreement with the observed data (Stairs et al. |2000| ), we model the changes of the 
pulse shape in a precession cycle with period ~ 1000 d and find that the variation of the 
pulse duration follows from a standard beam pattern in each cycle. 


Key words: pulsar: individual(PSR B1828-11) - stars: neutron 


1. Introduction 

Analysis of the long-term observation of the spin behavior of the isolated pulsar PSR 
B1828-11 reveals periodic variations both in the pulse shape and the slow-down rate of 
the pulsar and shows strong Fourier power at periods of ~ 1000,500,250, and 167 d 


(Stairs et al. 2000). Close correlations of the periodic changes in the pulse shape and 
duration with variations in the spin-down rate of the pulsar strongly suggest that the 
star’s spin axis is freely precessing around the star’s symmetry axis. 


Recently Link & Epstein ( |200lD studied the behavior of the observed pulse duration 
of PSR B1828-11 based on the free precession model. Since both 500 d and 250 d Fourier 
components have dominant contributions in the observed variations of period residual 
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Ap, its derivative Ap, and pulse shape of the beam, they suggested that the actual free 
precession period of the star is close to the strongest Fourier component ~ 500 d. Then 
a coupling of nearly orthogonal magnetic dipole moment to the star’s spin axis would 
provide the observed harmonic at period 250 d. On this basis, they modeled the pulse 
duration variations as a function of beam’s sweep angle A0 (see below for definition), 
and found that the beam pattern of radiations must be non-standard: the beam duration 
is larger for beam sweeps farthest from the dipole axis. As a result, both upper and lower 
parts of the emission region are wider than its middle, see Figure 2. As they mentioned 
this is required by the precession interpretation of PSR B1828-11. 

Recent theoretical investigations on the free precessing motion of PSR B1828-11 pro¬ 
vide new explanations for the reported data with the fundamental period 1000 d. 


which was originally suggested by Jones & Andersson ( 2001 ). Rezania (2002) consid¬ 
ered the case in which the magnetic field of the star varies with time while the star is 
precessing. Then he found a condition under which a coupling of the star’s crust with 
the time-varying magnetic radiation torque would produce the whole observed Fourier 
spectrum consistently. In a good agreement with the data he found that the fundamental 


precession period would be ~ 1000 d. Alternatively, Wasserman (|2002D showed that in 
general an oblique rotator must process. By analysis of the mechanical energy of the 
system, he found that the minimum energy state for such star is a state where the star 
processes. For strong magnetic stresses in the star’s type II superconductor core he es¬ 
timated the precession period as Ppre — 2460 d/(/3cosxi3i2iJi5), where B 12 = P/IO^^ 
G is the star’s magnetic field strength, iJis = iJ/10^® G corresponds to the first critical 
field strength in a type II superconductor {Her ~ 10^® G), /3 ~ I and x is the inclination 
angle of magnetic symmetry axis. For the case PSR B1828-11 with B 5 x 10^^ G, the 
precession period will be Ppre ~ 1000 d, if the inclination angle is y ~ 60°. 

In this paper we study the pulse shape of PSR B1828-11 by assuming that the actual 
precession period of the star is ~ 1000 d rather than ~ 500 d. We note that Link & 
Epstein’s calculations give the right behavior of the pulse shape variations provided the 
actual precession period is close to ~ 500 d only. For the case Ppre — 1000 d, their analysis 
gives an incorrect prediction for the observed shape variations. Here, we generalize link 
& Epstein’s analysis and show that, for Ppre — 1000 d, to get the correct behavior for 
the shape variations (compared with data) the radio beam of PSR B1828-11 must have 
the standard pattern. 


2. The model 

Let X, y, and z represent an inertial coordinate system 5, and an observer is in the x — z 
plane with a; > 0 and z > 0. Now consider a rigid, biaxial rotating star with angular 
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velocity u;, the principal axes 61 , 62 , 63 , and corresponding principal moment of inertia 
Ii = I 2 ^ Is- The star’s angular momentum L, along the z-axis in 5, is misaligned to the 
star’s symmetry axis 63 by a wobble angle 9, ie. L -63 = Lcos9, see Figure 1. To study the 
pulse shape variations we look at the variation of the polar angle of the beam with respect 
to L (fixed in S) that is equal to the polar angle of the magnetic dipole moment m.|^ The 
azimuthal and polar angles of the magnetic dipole moment m, d* and 0 in the inertial 
frame S are given by tan <i> = mylvrix and cos 0 = vrizlm = sin 9 sin ip sin y + cos 9 cos y, 
where 9 is the wobble angle, ip = tan“^ a;i/a ;2 = 7 r /2 — ojpt — /3 {oji = Ci ■ u>) and x 
is the inclination of magnetic field symmetry axis from the star’s symmetry axis 63 , ie. 
m ■ 63 = mcosx- Here ujp = 27r/Ppre is free precession frequency and /3 is a constant 
phase. Expanding cos0 to the first order in 9 one finds 0 = x “ 9cos{ujpt + /3). The 
latter shows the beam polar angle 0 changes sinusoidally about x, as the star precesses. 

Now let 7 = ^ + X be the polar angle of the observer in S (the constant angle between 
the observer and the angular momentum vector), see Figure 1. Following Link & Epstein 
( |2001 ) we define the sweep angle A0 as the difference in polar angle of observer, 7 , and 
the dipole, 0 at the time of the pulse: 

AQ = j — Q = ^ + 9cos{ujpt + f3) + 0(9'^). (1) 

Here ^ is a free parameter which will be fixed later by fitting the data. Link & Epstein 
assumed that the pulse duration ic is a function of A0 only, and has an extremum at 
A0 = 0. Then they expanded the pulse duration w in terms of A0 asw = wo + W 2 {AQ)'^. 
As it is clear, the latter expression would provide both Ppre and Ppre/2 components (due 
to the cos{u}pt + (3) and cos^{ujpt + (3) terms) in the pulse duration. As a result, it gives the 
correct behavior of the observed pulse duration for Pp^ — 500 d. But with Ppre — 1000 d, 
the 250 d Fourier component is missing. To get the correct behavior for the pulse shape 
variations during a 1000 d precession cycle, we take a more general expansion rather than 


Link & Epstein ( 2001 ) as 

re = rco + W2(A0)^ + ■u;4(A0)^. (2) 

It is interesting to note that equation (^ provides the contribution of Ppre, Tpre/2, 
Ppre/3, and Ppre/4 terms in the pulse duration. So with Ppre — 1000 d, one would expect 
to observe 500 d, 333 d, and 250 d Fourier components in the pulse shape parameter. 
As reported by Stairs et al. (2000 | ), the 1000 d and 333 d components though small 
(in comparing with 500 d and 250 d components) they have non-zero amplitude in the 
observed pulse shape parameter. These terms are missing in the Link & Epstein’s model. 
^ Note that though the pulsar beam’s direction is not necessarily in the same direction as the 


dipole moment m, for simplicity Link & Epstein (2001) defined a pulse as occurring when the 
azimuthal angle $ of the magnetic dipole equals to the azimuth of the observer. As a result, <i> 
is the observed pulse frequency. To compare our concluding result with one obtained by Link & 


Epstein (2001), we use the same definitions. 
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Without loss of generality we assume that = 1. Following Stairs et al. (2000), we 
define the shape parameter as S = where and Aw are the fitted heights of 

the narrower and wider standard profiles respectively, so that S cs 0 for the wide pulses 
and S' ~ 1 for narrow ones. As a result, one can relate the shape parameter S to the 
pulse duration w of the observed beam by 

w = max(rc)(l — S) -I- min(ri;)S, (3) 

where max(iy) and min(r(;) are the maximum and minimum values of the beam duration 
in the precession cycle. Combining equations (H) and (^, one finds for u >2 > 0 
r(;2 + (A0)2 / A0 


^ ^ «^2 + (|?|+ 0)2 
while for iy 2 < 0 we have 


l?l+0 


S/So = 


l«)2l-(Ae)^ 

i^2i-(i«i-r0= vi«i+»y 


Ae y 


if \W 2 \ > (ICI 


1 - 


l^2|-(mS)= (it+ej ^ 2 ! < (ICI 

where Sq is a normalization factor. It is clear that the shape parameter depends on cOp, 
f3, 9, and both the sign and magnitude of W 2 - With the given values of 6* ~ 3° and 
Ppre — 1000 d the shape parameter is determined by W 2 , and sign(r(; 2 ) completely. In 
a good agreement with data, we find that W 2 must be negative with magnitude larger 
but close to |?C 2 | > (1^1 + 9)'^- Inserting the latter results in equation (||), one can easily 
show that the pulse duration w is bigger when the beam sweeps closer to the dipole axis 
(A0 = 0) rather than when it sweeps farther (A0 ^ 0). This means the beam pattern is 
standard: the beam duration is smaller for beam sweeps farthest from the dipole axis. Our 


Ae y 


0)2, 

0)2, 


(4) 


(5) 


result is completely in contrast with the result of Link & Epstein (2001) who found that 
a non-standard beam pattern (the beam duration is larger for beam sweeps farthest from 
the dipole axis) is required to explain the data. In Figure 2 we sample the pulse profile 
for different observer’s viewing angles (the solid closed curve) and compare with the one 
(the dashed parabola) obtained by Link & Epstein ( ^001 ). We find the radiation pattern 
of PSR B1828-11 is very close to multipole field radiation pattern with {£,m) = (1,±1) 
where £ and m are orbital and azimuthal numbers. As the star precesses, the observer 
sees different sweeps corresponded to different viewing angles. Line B corresponds to the 
viewing angle (closest to the dipole axis) for which pulse is the widest. S' ~ 0 (narrowest, 
S ~ 1) provided the star precesses with precession period Ppre — 1000 d (Ppre — 500 
d). Accordingly, the pulse profile is narrower, S ~ 0.4 (wider, S ~ 0.4) for line A and 
narrowest, S ~ 1 (widest, S ~ 0) for line C. 

By choosing an appropriate values for \w 2 \ and ^ one can easily fit the data. We found 
an acceptable fit on the pulse shape data of PSR B1828-11 by choosing |■u; 2 | — 0.0034 and 
^ = —0°.01 for the given values of 0 = 3°.2 and ojp = 27r/1016 d“^, see Figure 3. Here 
the value of (|^| -I- 0)2 ~ 0.0031. We note that our fit curve which is calculated for free 
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precession period Ppre = 1016 d, is indistinguishable from the one that was calculated 


by Link & Epstein (2001) with free precession period Ppre = 511 d. 


3. Summary 


Previous study on the pulse shape variations of the isolated pulsar PSR 1828-11 by Link 
& Epstein ( ^001 ) is based on the free precessing star with precession period ~ 500 d. 
Recent theoretical studies (Rezania ^002 , Wasserman 2002) on the free precession of 
PSR B1828-11 suggested that the actual precession period might be ~ 1000 d as seen 
in observations. In this paper, by assuming that the actual free precession period of the 
isolated pulsar PSR B1828-11 is Ppre — 1000 d we studied the star’s pulse shape variations 
during a precession cycle. To get the correct behavior in the pulse shape variations in 
comparing with the observed data, we expanded the pulse duration w to forth order of 
A0 (the difference between polar angles of the observer and dipole axis at the time of the 


pulse). The forth order term is not considered by Link & Epstein ( 2001 ). We found that 
the pulse duration w is larger when the beam sweeps closer to the dipole axis (A0 = 0) 
rather than when it sweeps farther to the dipole axis (A0 ^ 0), see Figure 2. This means 
that the variation of the pulse duration follows from a standard beam pattern in each 


cycle. Our result is in contrast with the result of Link & Epstein (2001) who found that 
a non-standard beam pattern is required to explain the data. 

Further, our model provides the contribution other Fourier harmonics, 1000 d and 
333 d, in the pulse shape parameter as seen in data. Although these terms have smaller 
(but non-zero) amplitudes rather than 500 d and 250 d Fourier components, they are 
forbidden in the Link & Epstein’s model. In addition, by taking into account a term 
(A0)® in equation (||), one can explain the 167 d Fourier component seen in data of PSR 
B1828-11 (Stairs et al. pOOO ). The latter which corresponds to Ppre/6 Fourier component 
for Ppre — 1000 d is also forbidden for Ppre — 500 d free precession model. 

Finally, we note that new observations of the pulse shape variations of PSR B1828-11 


would be able to determine its radiation pattern. As predicted by Link & Epstein (|200l| ), 
in a 511 d precession cycle the observer would see both the broad upper and lower 
parts of the radio beam. This would be a critical observational evidence for the proposed 
theoretical models, since by assuming Ppre — 1000 d as the fundamental precession 
period, the middle of the beam will be broader. 
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Fig. 2. The beam pattern for the radio emission from PSR B1828-11. The dashed parabo¬ 


las represent the beam pattern proposed by Link & Epstein ( |2001[ ) with Ppre — 500 d, 
while the solid closed curve correspond to the 1000 d precession period. The radiation 
pattern is very close to multipole radiation pattern with {£,m) = (1,±1). The lines B, 
A, C represent the viewing angle at which the beam is the narrowest {S ~ 1), wider 
{S cs 0.4), and widest {S ~ 0), respectively, for Link & Epstein’s model, while corre¬ 
spond to the widest {S ~ 0), narrower {S cs 0.4), and narrowest {S ~ 1), respectively, 
for the model discussed here. 
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Fig. 3. The pulse shape parameter data S for PSR B1828-11 (from Stairs et al. 2000| ). 
The shape parameter is defined as S' = where and Aw are the fitted heights 

of the narrower and wider standard profiles respectively, so that S ~ 1 for the narrowest 
pulses and S ~ 0 for wider ones. The data points are obtained by averaging S over 
multiple bins, and the solid curve is the fit calculated with Ppre = 1016 d, see text. 





























